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Abstract 

The observed clustering of ultra-high energy cosmic rays suggests 
the existence of a neutral component. The models with violation of 
Lorentz invariance may explain this component by neutrons becom- 
ing stable above some threshold energy Eq. The protons, in turn, 
may become unstable above some energy E\ > Eq. We calculate the 
dependence of the threshold energies Eq and E\ on the parameters 
of the model and find E\/Eq J> 1.5. We argue that the characteristic 
threshold behavior of charged and neutral components may be used as 
the specific signature of models with violation of Lorentz invariance. 
The existence of the neutron stability threshold Eq can be investigated 
with already existing data. 

1 Introduction 

Ultra-high energy cosmic ray (UHECR) experiments are believed to be one of 
the promising places where physics beyond the Standard Model may be dis- 
covered. Until recently this expectation was based on a number of arguments 
of which the non-observation |l| of the GZK cutoff is most acknowledged. 
The case is strengthened by the absence of nearby candidate sources in the di- 
rections of the observed UHECR, large-scale isotropy of arrival directions and 
generic difficulties with acceleration of particles to energies of order 10 20 eV. 

The apparent absence of the GZK cutoff by itself can be explained within 
the conventional physics. For instance, several models have been proposed 
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which assume isotropization of arrival directions in strong magnetic fields 
and attribute the observed UHECR to neutron stars of our Galaxy , Virgo 
cluster 0] or the nearby radio-galaxy Cen A ||. Regardless of a particular 
model, assuming hard injection spectrum in combination with local over- 
density of sources at scales ~ 50 Mpc makes the discrepancy between the 
observed and predicted spectra significantly smaller [|, 0. 

A new ingredient in this problem is the small scale anisotropy of UHECR. 
The study of arrival directions of UHECR reveals the existence of clusters of 
events || || [K|. The analysis based on angular correlation function shows 
that the typical size of clusters is comparable with the experimental angular 



resolution [pi] , [12| . A most natural interpretation of this result is the exis- 
tence of point-like sources of UHECR whose image is smeared by errors in 
determination of UHECR arrival directions. 

If this interpretation is believed, it has several important consequences. 
First, it excludes the models which assume diffuse propagation of UHECR 
in strong magnetic fields. Second, it makes very unlikely the scenarios which 
explain the absence of the GZK cutoff by assuming hard spectrum and lo- 
cal overdensity of sources. Indeed, in the latter case the arguments based 
on statistics of clustering impose unrealistic bounds on the local density of 
sources [|13|, [14]]. Finally, the observed clustering suggests the existence of 



neutral primary particles, since otherwise they would be deflected in the 
Galactic magnetic fields and clusters would not be so tight. 

Even stronger argument in favor of neutral particles follows from the cor- 
relation of arrival directions with BL Lacertae objects (BL Lac) fll5]| . BL 
Lacs are located at cosmological distances. For instance, two of the BL Lacs 
which coincide with triplets of UHECR are both at ~ 600 Mpc. If primary 
particles of the triplet events indeed came from BL Lacs, they must be neu- 
tral; moreover, they must be able to propagate over cosmological distances 
without substantial attenuation. 

The existence of neutral particles imposes further constraints on possible 
models of UHECR. Within the Standard Model, there are two stable neutral 
particles, photon and neutrino. The photon attenuation length due to e + e~ 
pair production on infrared background is of order 10 — 20 Mpc in the energy 
range 10 19 — 10 20 eV [fL6|] , so photons cannot explain the observed events 
unless the infrared background is extremely low and/or their initial energies 
are extremely high. 

Neutrinos by themselves cannot be particles which initiate airshowers if 
standard neutrino cross sections are assumed. They have to be converted into 
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hadrons and photons through interactions with primordial neutrino back- 



ground 1 17, [18L [L9|. Photons obtained in this way could account for the 



observed tight clusters. There are generic problems with this scenario |]20 
First, the smallness of neutrino cross sections implies large neutrino flux. 
Second, the production of neutrinos by accelerated protons requires signifi- 
cantly higher proton energy, while even energies of order 10 20 eV are difficult 
to achieve in most acceleration sites. If these difficulties can be overcome, 
neutrino may be an appealing candidate. 

While the situation with neutrinos is unclear, it is worth considering other 
possibilities. Outside of the Standard Model, several candidates for neutral 
primary particles have been proposed. One of them is a light SUSY hadron 
2T|, which can be, e.g., a w<is-gluino bound state or gluon-gluino bound 



state. Another possible candidate suggested recently is a light sgoldstino 
The models of this type have an advantage of explaining the neutral 



component and the absence of the GZK cutoff at the same time. The disad- 
vantage is that primary particles have to be produced by accelerated protons, 
which leads to losses in both flux and energy. As a result, already very tight 
requirements on acceleration energy and efficiency become even tighter. 
An attractive model which shares the same advantages but does not have 



the latter drawback has been proposed in Ref. p4 |. It is based on possible 



violation of Lorentz invariance at high energies (see, e.g., |2B| for 

further discussion of Lorentz invariance violation in the context of UHECR). 
Models of this type can be constructed without giving up the main princi- 
ples of quantum field theory in the framework of the brane world scenario 



29[ . Under rather general assumptions, violation of Lorentz invari- 



ance can be described phenomenologically by introducing different maximum 



velocities for different particles. As noted in Ref. ||24|| , these parameters can 
be arranged in such a way that the main cause of the GZK cutoff, the pion 
photoproduction process, is not operative. At the same value of parameters, 
neutron may become heavier than proton at sufficiently high energy, and 
therefore stable. Thus, it can serve as a neutral primary particle which is 
indistinguishable from proton in other respects. Moreover, protons (which 
are stable at low energies) may be accelerated in a usual way. Once they 
reach the threshold energy and become unstable, they decay into neutrons 
by "/5-decay" with no loss in the flux and practically no loss in energy. An- 
other way to produce UHE neutrons is in collisions of UHE protons with 
synchrotron photons in the acceleration site. 

How the above models can be distinguished? A clear signature of the 
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neutrino model is the direct detection of ultra-high energy neutrinos, whose 
flux is expected to be in the range of sensitivity of the Pierre Auger and 
Telescope Array experiments . An indirect argument in favor of neutrino 
models would be identification of the neutral component as photons. On 
the contrary, the hadronic neutral component speaks in favor of light SUSY 
hadron models or violation of Lorentz invariance. 

The purpose of this paper is to point out the signature which can dis- 
criminate between the latter two possibilities. It is based on the fact that 
within the model of Ref. [Q the behavior of proton and neutron masses 
with energy follows a certain pattern: with increasing energy, first the neu- 
tron becomes stable at some energy E , and then proton becomes unstable 
at some higher energy E% > E . In other words, there is a window in which 
both neutron and proton are stable, while at higher energies proton becomes 
unstable. Thus, in the model of Ref. [M\ , one generically expects that there 
is only charged component at E < Eq, both charged and neutral components 
at E < E < Ei, and only neutral component at E > Ex. The appearance 
and disappearance of components is a threshold effect and is, therefore, a 
step-like function of energy. This step-like behavior is the specific signature 
of models with violation of Lorentz invariance. 



2 Kinematics of the nucleon decays 

Let us now calculate the width of the energy window in which the two com- 
ponents coexist. Following Ref. [21] consider the model characterized by 
the maximum attainable velocities of neutron, proton, electron and neutrino 
equal to c n , c p , c e and c u , respectively. Without loss of generality one may 
set 

c n = 1 • 

Certainly, c p , c e and c u should be very close to 1 as well. The differences 
between them are at the level of < 10~ 20 ||24| . As it was shown in Ref. [24|, 
the neutron is stable at high energies only if 

Cpi c e , c v ^ 1 . 

Moreover, if c p ^ 1 then proton becomes unstable at high enough energy. 
We will show that the ratio E /Ei is determined by the parameter 

c 2 -c 2 
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while the absolute magnitudes of both E and E\ are proportional to e \ 
where 

e^(cl-cl)^. (2) 

Following the formalism of Ref . J24J , let us find the energy E above which 
the neutron is stable, as a function of the parameters c p , c e and c u . In order to 
do this one has to find the minimum possible energy £' mn (p n ) of the potential 
decay products, i.e. proton, electron and neutrino, at fixed total momentum 
p n . The condition of stability of the neutron with the momentum p n and 
energy E n = v/p£ + m\ has the form 



VPl + m l < ^mn(Pn) • (3) 



As it was shown in Ref. p^| , the equality in Eq. (Q) may hold only at one 
value of the momentum p n . Consequently, the energy E is equal to 



E = V^o + < , (4) 
where p is a solution to the equation 



pl + ml = E mn (p ) . (5) 

Clearly, the minimum total energy of the proton, electron and neutrino with 
the fixed total momentum is achieved when the momenta of all three parti- 
cles are collinear. Thus, in order to calculate E mn {po) one should find the 
minimum value of the following function 



£n(p P ,Pe,Pu) = CpJpl + m 2 p c 2 p + c e ^/p 2 e + m\c\ + c v p u (6) 
of three positive variables p p ,p e and p v subject to the constraint 

P0= Pp+Pe+Pu ■ (7) 

This minimum may be either an extremum of the function S n (p p ,p e ,p u ), or 
belong to the boundary of the triangular region determined by the constraint 
(0) and positivity conditions 

Pp, Pe, Pu>0. 
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As is shown in the Appendix, the minimum is achieved at zero neutrino 
energy and momentum, 

p„ = 0. 

Therefore, in order to find the energy at which neutron becomes stable one 
should study the kinematics of two-body decay n — > pe. As a consequence, 
E does not depend on the neutrino velocity c u . The problem reduces to the 
following system of equations 

Po = Pe+ Pp (8) 
CePe CpPp 



m? p c 2 p 



Po + m l = Ce\/pl + m 2 e c 2 e + Cp^Jpl + m 2 p c 2 p . (10) 

Here eq.(§) is the momentum conservation condition (0) with p v = 0, eq.@ 
is the extremality condition for the function £(p p ,p e ,0) and eq.(^) is the 
stability condition @. It is straightforward to solve this system analytically 
in the with the result 



' m "" ( 2 mp + me)2 - 3.85 x 10 19 ( ^) eV, (11) 
c 2 p -l \ e 1 



in agreement with Ref. [24] . 

In the case c p ^ c e it is convenient to express the velocities in terms of 
the parameters a and e according to 

c 2 p = 1 + e 2 , c 2 = 1 + ae 2 , 

and rescale the momenta 

Qi 

Pt = -• 

e 

We assume e to be very small, but make no assumptions about the value of a. 



Expanding eqs.@-(10) to the lowest non-trivial order in e we obtain a system 
of polynomial equations for three unknowns q , q p , q e , which depends on one 
parameter a. This system can be solved for q numerically. The threshold 
energy is determined by the relation 

E = (12) 
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Figure 1: The dependence of the rescaled threshold momenta go ~ ^Eq and 
q\ ~ eE\ on the parameter a. 

where the difference between Eq and po has been neglected. The numerical 
solution for q as a function of a is plotted in Fig. [I]. 

Let us now find the threshold energy Ei above which the decay p — > nev 
is kinematicaly allowed. We follow the similar procedure as in the case of 
neutron, i.e., find the threshold momentum p\ from the equation 

c v\/pi + m l c l = E mp (pi) , (13) 

where E mp (p p ) is the minimum total energy of a neutron, electron and neu- 
trino with the fixed total momentum p p . This energy is obtained by the 
minimization of the function 

£ P (Pn,Pe,Pv) = VpI + ml + c e ^/p 2 e + m? e c 2 e + c u p v (14) 
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under the constraints 



Pn+Pe + Pu=Pp, Pn,Pe,Pu>0. (15) 

As in the case of a neutron one may show that the minimum is reached at 
p u = and therefore it is sufficient to study the kinematics of the two-body 
decay p —> ne when both momenta p p and p e are non-zero. The system of 
equations determining the threshold momentum p\ has the following form 
(cf. Eqs. (§)-©) 

Pi = Pe+Pn (16) 

r < P< /'» (17) 



p\ Pl + c l m l = c eVpl + C l m l + \ZpI + m l ■ ( 18 ) 



Again, the solution is straightforward to find in the case c p = c e (cf. Ref. @]), 



\/mi — (m r , — m e ) 2 10 / 10 12 \ 

E 1 ^^^ c ±- ~ 5.83 x 10 19 f —— \ eV, (19) 

At c p 7^ c e , the solution can be obtained numerically. The resulting function 
qi(a) is plotted in Fig. [I]. At small e, the threshold energy E 1 is related to 
q x {a) by 

E X = q -^. (20) 

Comparing Eqs.([Tl|) and ( |T9| ) one finds that the ratio Ei/E Q ~ 1.5 at 
c p = c e (i.e., at a = 1). It grows with a, as is shown in Fig. |^. The 
dependence becomes nearly linear at the values a ^ m p /m e ~ 2000. 



3 Discussion 

As we have seen in Sect. 2, the threshold energies of neutron stability and 
proton instability, Eq and E\, depend on two parameters e and a as defined 
by eqs. (O) and (]2"0|). The parameter e sets an overall scale of E and Ei, 



while the parameter a uniquely determines the ratio Ei/E . This ratio is 
always larger than ~ 1.5. 

In order to be of phenomeno logical interest, the neutron stability thresh- 
old E must be low enough. If clustering in Yakutsk data [B L) is attributed to 
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Figure 2: The ratio of threshold energies E\/Eq as the function of the pa- 
rameter a. 



neutrons, the existence of triplet with energies (2.5; 2.8; 3.4) x 10 19 eV implies 
E < 2.5 x 10 19 eV. At a <C m p /m e this requires in turn e > 1.5 x 1CT 12 , or 



c 2 - c 2 

p fi 



> 2 x 10 (this constraint is somewhat weaker for large values of 
a, see figfj). This constraint does not contradict the existing upper bounds 
|24|| on c 2 — 1 and c 2 — 1. Interestingly, the threshold at E may, in prin- 



ciple, be detected (or ruled out) with already existing experimental data at 
low energies E < 4 x 10 19 eV. For this analysis the key issue is the angular 
resolution of the experiment, so the AG AS A data seems the most suitable 
choice. Note that the improvement of upper limits on c 2 — c 2 may rule out 
the models which explain the neutral component of UHECR by the violation 
of the Lorentz invariance. 

The threshold of proton instability E\ is more difficult to detect since this 
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would require good identification of charged component. In principle, this 
can be done by making use of the deflections in the Galactic and extragalac- 
tic magnetic fields. Two different situations should be distinguished. If at 
energies smaller than E\ the extragalactic fields cause random deflections by 
angles larger than the typical distance between the sources, the charged com- 
ponent would look like a uniform background. It will be difficult, if possible 
at all, to separate such a component from unresolved sources. 

On the contrary, if random deflections are smaller than typical separation 
between the sources, one would see "halos" of the events formed by charged 
primary particles around clusters of the events with neutral primaries. If 



correlations with BL Lacs at large angles |15j are believed, current data 
favor the latter situation. If so, with enough statistics one will be able to 
identify the charged component reliably. Even better situation is possible if 
random deflections are negligible. In that case not only the charged particles 
can be identified, but their charges and the structure of magnetic fields can, 
in principle, be determined. 

It is worth noting that models with exotic neutral primaries also generi- 
cally predict the flux of charged particles, the direct protons from the source. 
The difference with the previous case is that in these models the proton com- 
ponent has a conventional GZK cutoff which occurs at fixed and calculable 
energy Eqzk, while the proton instability threshold Ei is, in general, differ- 
ent. If by chance Ei ~ Eqzk, the discrimination between the two models 
would be very difficult. 

To summarize, we argue that the characteristic threshold behavior of 
charged and neutral components in models with violation of Lorentz invari- 
ance can be used as their experimental signature. Upper bound on the neu- 
tron stability threshold E can be obtained from already existing data, the 
proton stability threshold E\ may be seen by the Pierre Auger experiment 
[31] where about 1000 events in the energy range 10" - 10 20 eV are expected 
in a year. 
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Appendix. Minimization of £ n (Pp,PeiPv) 

It is straightforward to check that function S n (p p ,p e ,p u ) has an extremum 
inside the triangular region if the inequality c p , c e > c u and the following 
condition is satisfied 

c u c e c u c p . . 

Pn > T~5 9M /9 me ' 7~9 9 \ 1 /9 m P ' V^-U 

[ci - ciyi 2 {c 2 p - ci) 1 ' 2 

The total energy of the proton, electron and neutrino in this extremum is 

Exipn) = c v p n + m p c p (c 2 p - c 2 u ) 1/2 + m e c e (c 2 e - c 2 u ) 1/2 . (22) 
Now one can readily check that the equation 

VpI + m l = Ei(Pn) 

is incompatible with the inequality ([Hj). Therefore, in order to determine the 
threshold momentum p Q of neutron stability one should find the minimum 
value of the function S n (p p , p e , p v ) on the boundary of momentum region. For 
that one should consider six cases when either one or two of the momenta 
p p , p e , p v are zero. 

Let us start with the case p e = 0. The function £ n (p p ,0,p v ) has an 
extremum with respect to p p and p u provided that c p > c v and 

The energy corresponding to this extremum is 

E 2 {p n ) = c u p n + m p c p (c 2 p - cl) l/2 + m e (? e . 

This energy becomes larger than the energy of neutron with the momentum 
p n at very low momenta p n > (m 2 l —m 2 )/2m e . Thus, this kinematical channel 
of the neutron decay is closed at energies higher than a few GeV even in the 
Lorentz-invariant case. It is straightforward to check that the same happens 
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in the = and when any two momenta of the decay products are 

zero. Therefore, the minimum of interest corresponds to the case p v = 0. 

The threshold of proton instability is found in a similar way The result is 
the following: in order to find the energy E\, where proton becomes unstable 
one should set p v = and study two-body decay p — > n, e. To illustrate 
some differences with the case of neutron decay, let us consider just one 
alternative kinematical regime, corresponding to point on the boundary of 
the region determined by ([15]) where only the proton momentum is non-zero. 
In this case equation fll3|) reduces to 

c P \Jp{ + m 2 p cl = ^Jpl + ml + m\c\ 

which gives 

2m e 



This value is indeed much larger than the one in the extremum with non-zero 
p e , see Eq. ( |T9"D and Fig. [TJ. 
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